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Aspects  of  Turbulent  Transport 

V.  Naulin  and  J.  Juul  Rasmussen 

Association  EURATOM-Ris0  National  Laboratory 

Turbulence  is  an  ubiquitous  phenomenon  in  plasmas,  which  have  usually  large  amounts  of 
free  energy  available  that  can  drive  the  turbulent  velocity  fluctuations.  Universal  properties 
of  turbulence  are  its  abilities  to  mix  and  to  transport.  The  turbulent  transport  is  often 
found  not  to  be  described  by  a transport  coefficient,  as  Fick’s  law  suggests.  We  discuss  the 
consequences  of  this  and  the  relationship  between  passive  particle  and  plasma  transport,  that 
is  mixing  and  transport. 


1.  Introduction 

A characteristic  feature  of  turbulence  is  the  ability 
to  disperse  and  mix  particles  and  heat.  This  process 
is  rather  complicated  even  under  idealized  conditions 
of  homogeneous  and  isotropic  turbulence  and  is  af- 
ter many  years  of  research  still  far  from  being  un- 
derstood in  general.  Covering  a variety  of  important 
areas  from  the  diffusion  of  pollutants  in  environmen- 
tal flows  to  the  particle  and  heat  diffusion  in  magne- 
tized plasmas  this  topic  connects  basic  research  with 
applications. 

For  magnetized  plasmas  it  is  agreed  that  the  en- 
hanced levels  of  cross-field  transport  (anomalous 
transport),  observed  in  a wide  variety  of  devices 
including  tokamaks  and  stellarators  is  due  to  low- 
frequency,  electrostatic,  micro-turbulence,  with  the 
E x B-drift  as  the  dominating  velocity.  A good 
candidate  to  understand  and  explain  the  anomalous 
transport  from  first  principles  is  drift-wave  turbu- 
lence [1,  2,  3].  We  investigate  turbulence  of  magne- 
tized plasmas  and  discuss  the  influence  of  structures 
on  the  transport. 

2.  Particle  Dispersion... 

The  mean  square  displacement  of  an  ensemble  of  par- 
ticles in  a velocity  field  u(x,  t)  is  given  by  [5]: 

(x2(t))  = 2 1 (u2)  J (i  - ^)  RL{r)dT  (1) 

with  ( u(x(t ))  • u{x{t'),t'))  = Rl{1  — t')  (u2)  . Here 
Rl  is  the  normalized  Lagrangian  velocity  correlation. 
For  very  short  times  the  correlation  is  close  to  unity 
in  the  so-called  ballistic  limit: 

= t2  (u2)  . (2) 

For  very  long  times  the  integral  evaluates  to  the  La- 
grangian integral  77 

t rOO 

Rl^Jt  ~ / RL{r)dT—.TL.  (3) 

Jo 

and  provided  that  77  is  finite  we  obtain 

(x2(t))  — ItTL  (u2)  = 2 tlD  (4) 

This  limit  ( t tl)  is  often  called  the  diffusion 

limit.  Note  that  tl  differs  from  the  Eulerian  correla- 
tion time  available  to  laboratory  measurements.  For 
homogeneous,  isotropic  turbulence  and  times  longer 


than  the  Lagrangian  correlation  time  the  particle  dif- 
fusion will  thus  be  normal.  However,  there  might  be 
large  intermediate  ranges  depending  on  the  correla- 
tion times  and  scales,  especially  if  the  correlations 
reach  the  spatial  or  temporal  extend  of  the  system. 

3.  ...  and  Transport 

The  question  if  test  particle  transport  is  connected 
to  the  turbulent  radial  transport  of  plasma 

r = (nvT) 

by  the  ExB  velocity  (vr  = — dy<f> ) is  far  from  trivial. 
We  consider  the  Hasegawa-Wakatani  equations 
(HWE)  [4]  as  minimal  model  for  drift-wave  type 
plasma  turbulence: 

dtn  + dy4>  + {<t>,  n}  = -C  ( n-<j>)  + /iV2n  , (5) 

dtV2  </>+{<!>, V2cj>}  = -C(n-<t>)  + pVi<t>  (6) 

The  deviation  from  adiabaticity,  given  by  the  param- 
eter 1/C,  leads  to  an  instability.  The  HWE  Eqs.  (5, 
6)  contain  two  limits:  i)  C — * 0,  where  the  density 
and  potential  decouples,  thus  the  system  describes 
standard  2D  Navier-Stokes  turbulence  with  n pas- 
sively advected.  ii)  C — ► 00,  which  corresponds  to 
the  adiabatic  density  response  (j>  ss  n and  the  HWE 
reduce  to  the  Hasegawa-Mima  equation  (HME). 

The  evolution  of  the  energy  E = ^ f ((V<p)2  + n2)d.V 
and  generalized  enstrophy  W ~ | f{a>  — n)2dV  is 
governed  by 

dtE  = J —ndyip—C  ( n — ^)2-g(w2+(Vn)2)d2r  ,(7) 

d/W  = J —ndy4>  — p(V(n  — u)))2d?r  . (8) 

Here  w = V x v = V2tp  is  the  vorticity.  The  driving 
term  is  the  integrated  E x B-flux 

To  = J Fd2r  = J nu<pd2r  = J ~ndyipd2r, 

which  mediates  the  instability  when  n (j>.  In  the 
inviscid  limit  the  HWE  has  a Lagrangian  conserved 
quantity  n = (w  — n + x).  We  employ  the  conserva- 
tion of  n to  estimate  the  radial  dispersion  of  the  fluid 
elements.  From  x — xo  = — (£  — £0) , where  / = uj  — n. 
We  obtain: 

({x  - x0 )2)p  = <£2)p  + (Co)P  - 2(£0£)p,  (9) 
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Figure  1:  Diffusion  coefficient  and  Flux  for  various 
values  of  C. 

where  the  average  (-)p  is  taken  over  all  fluid  elc- 
ments/particles.  If  the  turbulence  is  homogeneous, 
we  can  replace  the  averaging  over  the  fluid  elements 
with  averaging  over  the  fields  ((-)p  = (•)/).  In  the 
time  asymptotic  limit  we  may  then  assume  that  the 
last  term  on  the  rhs  vanishes  (this  designates  the 
correlation  between  the  fluid  PV  and  the  initial  fluid 
PV).  Taking  the  time  derivative  and  using  Eq.  (8) 
we  obtain: 

dt((x-xof)P~dt(e)f  = 2dtW.  (10) 

Finally  we  arrive  at 

dx  = r0 , (li) 

what  is  confirmed  by  numerical  experiment  see  Fig- 
ure 1. 

4.  Transport  PDF 

In  models  where  the  fluctuations  are  driven  by  a local 
instability  and  the  density  fluctuation-level  is  small 
compared  to  the  background  density  the  probability 
distribution  function  (PDF)  of  as  well  density  and 
radial  velocity  is  given  by  a normal  distribution. 
When  evaluating  for  the  PDF  of  T we  note  that  the 
folding  of  two  normal  distributions  is  not  a normal 
distribution.  Moreover  if  n and  <j>  are  correlated  the 
average  of  the  PDF  may  become  non-zero  and  lead 
to  a net  transport. 

If  the  transport  PDF  is  evaluated  not  for  the  lo- 
cal values  of  the  transport,  but  averaged  over  some 
domain,  the  transport  PDF  has  to  converge  to  a 
stable  distribution.  If  the  second  order  moment  is 
finite,  the  central  limit  theorem  tells  us  that  this 
is  a normal  distribution.  If  the  second  order  mo- 
ment does  not  exist,  than  we  have  a Levy-Parcto 
stable  PDF  and  thus  expect  a power  law  tail  in  the 
PDF,  as  long  as  all  the  averaging  procedures  are  car- 
ried out  for  times/lengths  well  below  the  correlation 
length/time-scales.  A power  law  tail  in  the  PDF  of 


the  averaged  flux  reflects  thus  the  existence  of  long 
range  correlations  in  the  system,  respectively  that 
the  turbulence  is  intermittent.  Correlations  that  can 
manifest  themselves  via  the  occurrence  of  coherent 
structures  within  the  turbulence. 

5.  Turbulence  and  Structures:  The 

quest  for  correlations 

Two  dimensional  turbulence,  as  the  turbulence  of 
magnetized  plasmas,  is  characterized  by  a dual  cas- 
cade: of  enstrophy  to  the  small  scales  and  energy  to 
large  scales.  Moreover  if  the  vorticity  is  a function  of 
the  potential,  then  the  convective  non-linearity  van- 
ishes and  thus  turbulent  de-correlation  mechanisms 
are  no  longer  very  effective.  The  regions  of  reduced 
nonlinearity  and  correlated  potential  and  vorticity 
are  also  known  as  coherent  structures.  These  struc- 
tures form  out  of  the  turbulence  and  can  survive  for 
long  times  compared  to  the  turbulence  de-correlation 
time.  Tracking  particles  in  the  turbulence  and  stat- 
ing if  they  arc  trapped  inside  structures  we  could  find 
that  this  introduces  an  intermediate  regime  for  the 
transport  [6]. 

6.  Turbulent  Equipartition 

Usually  turbulence  takes  place  in  an  confined  vessel 
with  boundary  conditions  and  the  plasma  is  contin- 
uously driven.  The  arguments  used  above  are  then 
no  longer  valid.  Indeed  based  on  the  conservation  of 
Lagrangian  invariants  we  could  show,  that  the  turbu- 
lence sets  up  characteristic  gradients,  independently 
of  the  magnitude  of  the  flux  generated.  Clearly  this 
breaks  Fick’s  law.  Moreover  it  leads  to  a bursty  type 
of  transport,  where  plasma  is  pushed  outwards  in 
large  structures  [7]. 

7.  Summary 

It  is  found  that  as  long  as  the  spatial/temporal 
extend  of  the  plasma  is  close  to  correlation 
times/lcngth  scales  the  turbulent  structures  induce 
a power  law  tail  in  the  PDF  of  the  transport.  This 
transport  is  then  not  described  by  Fick’s  law  on  dif- 
fusion and  turbulent  profiles  might  be  sustained  at 
various  levels  of  flux  through  the  system. 
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